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ABSTRACT 


This  report  sets  up  a  model  for  the  intensity  of  infrared  radiation  emitted  by  a  particulate 
exhaust  plume  generated  by  a  pyrotechnic  flare  burning  in  a  vacuum.  The  model  assumes  that  the 
exhaust  plume  is  optically  thin,  that  it  is  composed  mostly  of  discrete  particles,  that  the  particles 
have  a  size  spectrum  described  by  a  log-normal  probability  density  function,  that  they  cool  off 
entirely  by  the  emission  of  thermal  radiation,  and  that  they  all  travel  at  the  same  average  velocity. 
The  model  is  used  to  predict  how  a  magnesium-Teflon  exhaust  plume  would  look  when  viewed  as  an 
approximate  point  source  by  a  distant  infrared  sensor  and  also  to  analyze  the  data  acquired  from 
three  separate  magnesium-Teflon  flares  burned  in  a  large  vacuum  chamber. 
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A  MODEL  FOR  THE  INFRARED  RADIANCE  OF 
OPTICALLY  THIN,  PARTICULATE  EXHAUST  PLUMES 
GENERATED  BY  PYROTECHNIC  FLARES 
BURNING  IN  A  VACUUM 


1.  INTRODUCTION 


For  the  last  two  years  the  Air  Force  has  sponsored  work  in  Group  95  at  Lincoln  Laboratory 
aimed  at  understanding  the  characteristics  of  infrared  (IR)  radiation  emitted  from  the  exhaust 
plumes  of  pyrotechnic  flares  burning  in  a  vacuum.  Over  this  period  of  time,  a  model  has  been 
developed  which  has  proven  to  be  reasonably  successful  at  explaining  the  results  of  laboratory 
experiments.  This  model  assumes  that  the  vacuum  exhaust  plumes  are  optically  thin,  that  they 
are  composed  mostly  of  particles,  and  that  the  particles  cool  off  entirely  by  the  emission  of  ther¬ 
mal  radiation  as  they  travel  down  the  plume.  The  particles  are  assumed  to  come  in  a  variety  of 
sizes,  with  a  size  distribution  which  obeys  a  log-normal  probabiUty  density  function,  and  to  travel 
at  the  same  average  velocity.  The  model  will  work  for  any  exhaust  plume  composed  primarily 
of  discrete  particles,  but  it  has  been  applied  here  only  to  Mg-Teflon  flares,  because  working  with 
this  type  of  pyrotechnic  proved  to  be  experimentally  convenient. 

Following  this  introduction.  Section  2  outlines  the  mathematical  development  of  the  model 
equations  and  describes  the  basic  approximations  that  went  into  their  derivation.  In  particular,  a 
phenomenological  approximation  (instead  of  Mie  theory)  is  used  to  describe  the  variation  of 
exhaust  particle  grey-body  emissivity  with  particle  size,  aUowing  the  model  to  acknowledge  both 
the  geometric  irregularity  and  the  lack  of  chemical  purity  of  the  Mg-Teflon  exhaust  plume  parti¬ 
cles.  In  Section  3  the  model  is  simplified  further  by  assuming  that  all  the  exhaust  particles  have 
the  saine  size  ^  well  as  the  same  average  velocity.  The  resulting  savings  in  computer  time  made 
it  practical  to  investigate  how  the  exhaust  plume  of  a  Mg-Teflon  flare  would  look  when  viewed 
as  an  approximate  point  source  by  a  distant  IR  sensor.  Since  the  exhaust  plume  particles  were 
all  the  same  size,  and  this  size  was  adjustable,  it  was  also  a  straightforward  matter  to  investigate 
the  effect  of  parricle  size  on  plume  radiance.  Section  4  again  represents  the  particles  by  a  log¬ 
normal  distribution,  and  the  predictions  of  the  model  are  compared  to  IR  data  acquired  from 
three  separate  burning  flares.  These  flares  were  part  of  a  longer  series  of  tests  conducted  in 
March  1988,  which  were  designed  to  measure  the  radiance  of  Mg-Teflon  flares  burning  in  a 
vacuum.  The  concluding  section  summarizes  the  results  of  the  report. 
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2.  DERIVATION  OF  THE  MODEL  EQUATIONS 


When  measuring  the  size  distribution  of  small  irregular  particles,  it  is  customary  to  assume 
the  particles  are  spheres  and  then  allow  for  their  irregularity  by  introducing  undetermined  con¬ 
stants  into  the  formulas  for  the  surface  area  and  volume  of  a  sphere  [1],  Suppose,  for  example, 
the  size  distribution  of  a  collection  of  microscopic  cubes  with  sides  of  lengths,  were  examined  by 
some  technique  which  measured  directly  the  surface  area  of  each  particle.  Not  knowing  that  the 
particles  were  cubes,  each  particle  would  then  be  assigned  an  equivalent  spherical  radius  R  using 
the  rule  that  47rR2,  the  surface  area  of  a  sphere,  must  equal  the  surface  area  of  a  cube.  It 
•follows  that  j,  =  R(27r/3)*/2,  but  the  volume  of  that  same  particle,  if  it  were  truly  a  sphere,  would 
be  (4/3)7rR3.  Since  it  is,  in  fact,  a  cube,  its  volume  in  reality  is.i.^  =  (27r/3)^/2R3  =  fy(4/3);rR3 
where  fy=  (7r/6)’/2.  it  is  not  surprising  that  fy  <1,  because  a  sphere  can  be  defined  as  that  three- 
dimensional  closed  surface  which,  for  a  fixed  amount  of  area,  encloses  the  maximum  possible 
volume.  Therefore,  other  surfaces  can  be  expected  to  enclose  smaller  volumes. 

The  exhaust  particles  generated  by  the  Mg-Teflon  flares  described  in  this  report  have  had 
their  size  distributions  measured  by  a  modified  Malvern  model  2600c  particle  sizer.  The  Malvern 
particle  sizer  measures  the  amount  of  light  forward-scattered  off  microscopic  particles  illuminated 
by  a  helium-neon  laser.  From  this  it  can  estimate  the  fraction  of  particles  at  any  given  size, 
where  the  size  of  the  particle  is  taken  to  be,  in  the  case  of  spheres,  the  particle  diameter.  For  a 
collection  of  irregular  particles,  estimates  of  both  surface  area  and  volume  go  into  the  Malvern 
calculation  of  the  particle  size  distribution.  Therefore,  what  we  have  done  is  to  divide  the 
Malvern  size  estimate  —  that  is,  the  effective  spherical  diameter  of  a  particle  —  by  two  and  call 
the  resulting  number  the  particle  scattering  radius  R.  Then  the  surface  area  S  and  volume  V  of  a 
particle  of  scattering  radius  R  are  defined  to  be 

S  =  fs-47rR2  (2.1a) 

V  =  fy-y7rR3  (2.1b) 

where  fj  and  fy  are  unknown  positive  constants  independent  of  R.  In  effect,  all  the  particles 
inside  the  exhaust  plume  have  been  assigned  the  same  average  irregular  shape.  The  dimensionless 
constants  f^  and  fy  must  be  0(1)  in  magnitude  because  the  exhaust  particles  of  the  flare  have 
been  caught  on  glass  slides,  examined  under  a  microscope,  and  found  not  to  be  outrageously 
nonspherical  in  shape. 

In  the  model  developed  below,  it  is  assumed  that  the  exhaust  particles  of  the  burning  flare 
can  be  described  by  a  particle  size  distribution  f(R),  where  f(R)dR  is  the  probability  that  a  parti¬ 
cle  picked  at  random  from  the  flare  exhaust  plume  has  a  scattering  radius  between  R  and 
R  +  dR.  It  is  also  assumed  that  f(R)  is  a  log-normal  distribution 
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where  o  and  Ria  are  positive  constants  whose  values  will  be  measured  by  the  Malvern  particle 
sizer. 


In  order  to  predict  plume  radiant  intensity  away  from  the  flare,  it  is  necessary  to  know  how 
the  exhaust  particles  cool  off  as  they  travel  down  the  plume.  The  particles  in  the  exhaust  plume 
are  irregularly  shaped  chunks,  so  there  is  no  practical  way  to  calculate  their  emissivity  using  elec¬ 
tromagnetic  theory.  It  is  expected,  however,  that  when  the  particles  are  much  larger  than  the 
radiation  wavelength  X,  their  surface  emissivity  will  approach  a  constant  value  «  regardless  of 
their  shape;  and  that  when  they  are  much  smaller  than  the  wavelength,  that  is,  in  the  Rayleigh 
limit,  their  surface  emissivity  will  be  proportional  to  R/X  =  Kvl c,  where  R  is  the  particle  scatter¬ 
ing  radius,  v  the  radiation  frequency,  and  c  the  velocity  of  light  [2].  Therefore,  it  is  assumed  that 
the  exhaust  particles  have  an  emissivity  function  ^  which  obeys  the  rule 


for  R  vjc  >  l/27r 


(2.3) 


2;rR€j' 


c 


for  0^  Rv/c^  ll2n 


The  emissivity  of  a  microscopic  particle  which  is  both  chemically  pure  and  shaped  like  a 
sphere,  so  that  it  has  both  a  well-defined  complex  index  of  refraction  n  and  a  mathematically 
tractable  geometry,  can  be  predicted  using  electromagnetic  theory  (often  called,  in  this  context, 
Mie  theory).  Figures  2-1  and  2-2  graph  the  Mie-theory  emissivity  of  spheres  for  a  variety  of  dif¬ 
ferent  n  values.  In  each  case,  the  emissivity  changes  from  being  (more  or  less)  constant  to  being 
proportional  to  R/X  at  a  value  of  R  «  X/27r  =  dlnv.  Although  the  flare  exhaust  particles  are,  in 
general,  neither  chemically  pure  nor  shaped  like  spheres,  this  behavior  of  the  Mie-theory  emissiv¬ 
ity  suggests  that  K  =  cl  2irv  [see  Equation  (2.3)]  is  a  plausible  estimate  of  where  the  particle  emis¬ 
sivity  changes  over  from  Rayleigh  behavior  (proportional  to  R/X)  to  grey-body  behavior  (equal 


to  a  dimensionless  constant  c). 


Assuming  that  the  flare  exhaust  particles  are  isothermal  and  cool  off  entirely  by  the  emission 
of  radiation,  Ijot  is  defined  as  the  amount  of  radiant  power  per  unit  surface  area  emitted  by  a 
particle  of  temperature  T  and  scattering  radius  R.  It  follows  from  Planck’s  black-body  radiation 
law  and  Equation  (2.3)  that 


^tot  “ 


27rRe 


et,(T)d. 


(2.4) 


where 


eb.(T) 


c 

h 

k 


2irhv^ 

c2(ehj'/kT  _  1) 

speed  of  light  =  2.998  X  10®  m/s, 

Planck  constant  =  6.625  X  lO'^'^  erg-s,  and 
Boltzmann  constant  =  1.380  X  10**^  erg/®K. 
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EMISSIVITY 


Figure  2-1  Graph  of  Mie- theory  emissivity  versus  27t  X  (particle  radius)! 
(radiation  wavelength)  for  several  values  of  the  complex  index  of  refraction  n. 


IN  ALL  CASES 


X 


Figure  2-2  Graph  of  Mie-theory  emissivity  versus  27r  x  (particle  radius)! 
(radiation  wavelength)  for  several  values  of  the  complex  index  of  refraction  n. 


The  function  is,  of  course,  the  amount  of  radiant  power  emitted  from  the  surface  of  a 

black  body  at  temperature  T  per  unit  surface  area  per  unit  frequency  interval  at  frequency  v. 
Define  a  new  variable  ^  =  hi'/kT  and  rewrite  Equation  (2.4)  as 


Itot  =  T5 


/47r2RekA  ph 

\  hV  7*^0 


P^cad^  (2.5) 


hc/27rRkT 


where 

P/(a  =  ^^/(e^-l)  • 

The  integrals  in  Equation  (2.5)  cannot  be  simplified  unless  P/a  is  replaced  by  an  approxi¬ 
mation  P^a(a  such  that 

P/a(a  =  f^e-f  +  ^2e-^^  . 

When  n  >  1,  both  P^  and  P/a  have  the  limiting  form  for  ^  -  0  we  have 

«  1  +  X  +  0(x2)  and  thus  both  P^  and  P^  have  the  Umiting  form  a  +  ©(a)  when  f  -  0.  Ihe 
value  of  n  >  1  in  Equation  (2.6)  can  be  chosen  to  give  the  best  least-squares  fit  of  P^  to  P^. 
When  this  is  done,  SI  =  1.586,  and  Figure  2-3  shows  that  these  two  functions  then  match  each 
other  very  closely  indeed.  Figure  24  shows  that  the  maximum  error  involved  in  this  approxima¬ 
tion  is  several  percent  at  most,  which  is  quite  acceptable  for  the  purposes  of  this  report. 


Figure  2-3.  Graph  of  and  its  approximation  P/a  versus  f 
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Figure  2~4.  Graph  of  the  absolute  error  between  and  Piq  versus  f. 


Replacing  in  Equation  (2.5)  with  in  Equation  (2.6)  and  integrating  gives 
1  .^./27r€kA  _i4n-RkT\ 


where 


^(x)  =  24x(l-e-l/*)  +  6xn-4(I-e-«/x) 
-e-l/’^  (1/x2  +  6/x+  18) 
_ft-2e-n/x  (i/x  +  4/n) 

with 


(2.7) 


ft  =  1.586 


An  isothermal  particle  of  density  p  and  constant  heat  capacity  Cp  will  lose,  when  its  temper¬ 
ature  drops  by  an  amount  AT,  an  amount  of  thermal  energy  AE  given  by 


AE  =  fv  • 


4 

3 


ttR^p  Cp  AT 


(2.8) 


where  the  factor  fy  comes  from  the  formula  for  the  particle  volume  given  in  Equation  (2.1b). 
According  to  Equations  (2.7)  and  (2.1a),  the  particle  will  radiate  energy  at  the  rate  of 


dE 

—  =  -fs-47rR2ltot  =  -fsR2T4 


(^)i 


27rRkT) 


he 


(2.9) 
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Substitute  Equation  (2.8)  into  (2.9)  to  get 

l2Trm\ 

RpCph3c2 


dt 


f 


he 


(2.10) 


where 


\  =  emissivity  shape  factor  =  efg/fv 

Because  e,  and  fy  are  all  0(1)  quantities,  it  follows  that  is  also  0(1). 

In  Figure  2-5,  r  is  defined  as  the  distance  between  the  flare  and  a  particle  inside  its  exhaust 
plume.  Assuming  that  all  the  particles  inside  the  exhaust  plume  travel  (in  the  vacuum)  at  a  con¬ 
stant  velocity  v^.  Equation  (2.10)  can  be  rewritten  as 


dT 

dr 


67rAek4T4  ^  /27rRkT^ 


RpVoCph3c2 


he 


Equation  (2.11)  has  the  boundary  condition 


r  =  0 


=  T„ 


(2.11) 


(2.12) 


Figure  2-5  Steady-state  exhaust  plume  generated  by  a  flare  burning  in  a 
vacuum. 
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where  is  the  temperature  at  which  all  the  exhaust  plume  particles  leave  the  burning  flare. 
Equations  (2.11)  and  (2.12)  taken  together  define  a  function  T  such  that  T(r,R)  is  the  temp)era- 
ture  of  an  exhaust  plume  particle  of  scattering  radius  R  located  a  distance  r  from  the  burning 
flare. 

The  best  way  to  examine  the  behavior  of  function  T(r,R)  is  to  nondimensionalize  Equa¬ 
tions  (2.11)  and  (2.12). 


where 


dx  . 

—  =  -xi^(x) 
dq 


X  = 


^27rRk\ 

V  he  / 


(2.13a) 


(2.13b) 


/__3AgCk__\ 

\47r2R4p  v^Cp  y 


and 


(2.13c) 


(2.13d) 

Figures  2-6  and  2-7  give  x  as  a  function  of  q  for  two  different  ranges  of  Xq  values. 

The  steady-state  combustion  rate  of  a  burning  flare  can  be  represented  by  a  constant  P  such 
that  PAt  is  the  amount  of  mass  burned  by  the  flare  in  time  At.  Figure  2-5  shows  that  the  mass 
of  all  the  exhaust  particles  emitted  by  the  flare  in  time  At  must  lie  within  a  distance  v^  At  of  the 
flare’s  burning  surface.  Hence, 


q=0 


=  x„ 


27rRkT„ 


he 


q 

Figure  2-6  Graph  of  x  versus  q  for  =  0.25,  0.50,  0.75,  1. 
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Figure  2-7  Graph  of  x  versus  q  for  Xq  =  L  2,  5,  8, 

PAt  =  (Ap  VoAt)  •  •  ( j  TT  R3p)  f(R)  dR  (2.14) 

where  Ap  is  a  positive  constant  such  that  Apdr  is  the  number  of  particles  located  between  r  and 
r  +  dr  inside  the  flare  exhaust  plume,  p  is  the  average  density  of  the  exhaust  plume  particles  [see 
Equation  (2.8)],  fy  has  been  defined  in  Equation  (2.1b),  and  f(R)  has  been  defined  in 
Equation  (2.2). 

Equation  (2. 14)  can  be  rearranged  to  give 


A  3P 

P  47rpfv  Vo 


R3f(R)  dR  -•  = 


3Pe 


,-9a2/2 


4n-pfy  Vo  RI 


(2.15) 


*v  »o  ^la 

where  in  the  second  step  of  Equation  (2.15)  we  have  substituted  from  Equation  (2.2)  to  get 

R3  f(R)  dR  =  R^  .  (2. 16) 


A  radiometric  measurement  is  used  (see  Figure  2-5)  to  capture  the  radiation  emitted  by  all 
the  exhaust  plume  particles  located  between  r  and  r  +  dr  inside  the  exhaust  plume  of  the  burning 
flare.  Although  a  single  particle  will  have  a  distinct  irregular  shape  and  thus  cannot  emit  radia¬ 
tion  isotropically,  it  is  safe  to  assume  that  a  large  collection  of  such  particles  will,  being  ran¬ 
domly  aligned  with  respect  to  the  radiometer,  emit  their  collective  radiation  isotropically.  Thus, 
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on  the  average,  the  amount  of  radiant  power  per  unit  wavelength  at  wavelength  k  emitted  by  a 
single  particle  of  scattering  radius  R  and  temperature  T  into  the  solid  angle  Aftj  subtended  by 
the  Inframetrics  radiometer  is 

Afls 

- (f,47rR2)-^(R/X).ebx(D 

hTT 

where 


C] 

ebx(T)=  —7 - 

X5(e<=2/AT_  ]) 

Cj  =  3.74  X  10**^  w-m2  ,  and 
C2  =  1.44  X  10-2  ni-®K 


Function  is  defined  in  Equation  (2.3)  and  fj  is  defined  in  Equation  (2.1a).  Function  e5x(T)  is 
the  amount  of  radiant  power  emitted  from  the  surface  of  a  black  body  at  temperature  T  per  unit 
surface  area  per  unit  wavelength  interval  at  wavelength  X.  Call  G(X,r)  dr  dX  dfl  the  amount  of 
radiant  power  emitted  by  all  particles  located  between  r  and  r  +  dr  at  wavelengths  between  X  and 
X  +  dX  into  a  solid  angle  dfl.  If  the  exhaust  plume  is  optically  thin*  so  that  the  particles  do  not 
absorb  or  scatter  significant  amounts  of  radiation,  it  follows  that 


G(X,r)  =  f3  •  Ap  •  r®®  dR  R2f(R)  etx  (T(r,R))  -^(R/X) 
0 

^  3AePCi  e-9^^/2  mUf(R/k) 

47rpVo  R2^  X5  0  e^2/'*^T(r,R)  _  1 


(2.17) 


where  in  the  second  step  of  Equation  (2.17)  we  have  substituted  the  expression  for  Ap  in  Equa¬ 
tion  (2.15)  and  in  Equation  (2.10).  We  have  also  taken  the  constant  c  outside  the  integral  by 
defining  that 


C^(X) 


1 

27rx 


for  X  >  l/27r 
for  0  < x^  I/Itt 


(2.18) 


Equations  (2.2),  (2.11),  (2.12),  (2.17),  and  (2.18)  define  the  model  used  to  predict  the  radiant 
intensities  of  exhaust  plumes  generated  by  Mg-Teflon  flares  burning  in  a  vacuum. 


*  A  TV  camera  was  mounted  inside  the  space  chamber  to  watch  both  the  flare  and  its  associated 
exhaust  plume.  The  burning  surface  of  the  flare  could  easily  be  seen  even  when  viewed  through 
the  central  portions  of  the  plume,  showing  the  plume  to  be  optically  thin  to  visible  light.  IR 
radiation,  having  a  longer  wavelength,  should  be  attenuated  even  less. 
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3.  INVESTIGATION  OF  AN  EXOATMOSPHERIC  FLARE 


It  is  instructive  to  simplify  the  exhaust  plume  model  by  taking  the  limit  as  o  —  0  in  Equa¬ 
tion  (2.2),  transforming  f(R)  into  the  dirac-delta  function  6(R-Rja).  This  gives  all  the  exhaust 
plume  particles  the  same  scattering  radius  R^  =  Rja  and  allows  investigation  (without  using  an 
excessive  amount  of  computer  time)  of  how  the  exhaust  plume  of  a  vacuum-burning  Mg-Teflon 
flare  looks  when  viewed  as  an  approximate  point  source  by  a  distant  IR  sensor. 

Assume  that  the  particle  concentration  C(r,0)  (with  units  of  particles  per  unit  volume)  at  a 
given  point  inside  the  exhaust  plume  is  a  function  of  both  r  and  6,  the  spherical  polar  coordi¬ 
nates  of  that  point  as  defined  in  Figure  3-1.  Say  that 


C(r,0) 


forO^e^^n^x 
for  e  > 


(3.1) 


where  is  a  constant  having  units  of  inverse  length.  The  constant  Ap  has  already  been  defined 
by  the  requirement  that  Apdr  be  the  number  of  particles  located  between  r  and  r  +  dr  inside  the 
flare  exhaust  plume  [see  Equation  (2.14)  et.  seq.].  Test  whether  Equation  (3.1)  has  left  Ap  a  con¬ 
stant  by  writing 


max 


dO  sin  6  C(t,6) 


(3.2) 


Figure  3-1.  Three  coordinate  systems  used  to  investigate  how  a  Mg-Teflon  exhaust  plume 
would  look  when  viewed  by  a  distant  IR  sensor. 


13 


Substitute  Equation  (3.1)  into  (3.2)  and  get 
Ap  =  ^Afl 

where 


Afl  =  solid  angle  subtended  by  the  exhaust  plume  cross  section 
+  27r(l  -  COS0n,ax) 


(3.3) 


Substitute  Equation  (3.3)  into  (2.15),  take  the  limit  as  a  —  0,  and  replace  Rjj,  by  R,,  to  get 

^  =  A_/An=  - -  .  (3.4) 

47rpfv  Vo  r3  An 

Equations  (3.3)  and  (3.4)  demonstrate  that  Ap  is  indeed  still  a  constant,  so  Equations  (2.17)  and 
(2.18)  still  hold.  Take  the  limit  of  Equation  (2.17)  as  <t  —  0  and  set  Rj,  =  to  get 

3A,Ra«R<,/X)  /  C,  \ 

G  (X,  r)  = - —  — - )  .  (3.5) 

47rpVoRoX5  _  i  j 

The  subscript  “o”  has  been  added  to  G  as  an  indication  that  Equation  (3.5)  should  only  be  used 
when  all  exhaust  plume  particles  have  the  same  scattering  radius  Rq. 

Since  the  particles  traveling  down  the  flare  exhaust  plume  will  never  cool  down  to  absolute 
zero  in  a  finite  amount  of  time,  the  exhaust  plume  can,  in  principle,  be  regarded  as  infinitely 
long.  To  treat  the  plume  as  an  approximate  point  source  when  viewed  by  a  distant  IR  sensor,  we 
must  give  the  plume  some  finite  length.  This  is  done  by  picking  a  temperature  Tf<To,  where  Tq 
is  the  temperature  at  which  the  particles  leave  the  burning  flare,  such  that  all  particles  with 
temperatures  T  ^  Tf  are  assumed  to  contribute  only  negligible  amounts  of  radiation  to  the  total 
output  of  the  exhaust  plume.  (In  the  work  done  below  Tf  =  200®K.)  The  length  of  the  exhaust 
plume  corresponding  to  Tf  is  called  Lf  where  from  Equation  (2.11) 


/i)RoVoCph3c2  dT 

- ^ -  rio - 

67rAek4  ‘'t,  T4^(27rRokT/hc) 


with  R  =  Rq  because  now  all  the  exhaust  plume  particles  have  the  same  scattering  radius.  In  fact, 
Lf  is  the  maximum  value  of  the  radius  coordinate  r,  and  for  any  temperature  value  T  between  Tf 
and  To  we  can  associate  a  unique  radius  value  r  such  that 


r(T)  = 


pRoVoCph3c2 

6rrAjk^ 


dT' 


T  T'4^(27rRokT'/hc) 


(3.7) 


Approximate  in  Equation  (3.7)  by  the  function 


^a(x)  = 


X  +  wx2 
U_+ V(x  +  wx2) 


(3.8) 
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The  values  of  u  and  v  come  from  the  requirement  that 


^(x) 


X«1 


=  24+  —s  24.95, 


so 


and 


u  =  (24.95)-!  =  0.04008; 


=  6+  — s  6.501, 

x»l  fl3 


SO 

v  =  (6.501)-!  =  0.1538. 

For  the  values  of  and  T  typical  of  vacuum-burning  Mg-Teflon  flares,  it  is  known  that 
X  =  27rR(,kT/hc,  the  argument  of,^,  lies  somewhere  between  0.005  and  10.  Therefore,  in  Equa¬ 
tion  (3.8)  w  is  adjusted  to  give  the  best  least-squares  fit  of  to  ^  for  0.005  <  x  ^  10,  and 
w  =  44.02.  Figure  3-2  demonstrates  how  closely matches,^ for  x  values  between  0.005  and  10, 
and  Figure  3-3  graphs  the  absolute  error  between  and  as  a  function  of  x.  The  maximum 
error  never  exceeds  10  percent,  which  is  adequate  for  our  purposes. 


Figure  3-2  Graph  of  ^  and  its  approximation  versus  x. 
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Equation  (3.10)  is,  computationally,  a  much  quicker  formula  for  r(T)  than  is  Equation  (3.7),  and 
that  is  why  function  has  been  used  to  approximate 

To  calculate  the  response  of  the  IR  sensor  to  the  radiating  exhaust  plume  of  the  flare, 
simply  add  the  radiation  emitted  by  all  particles  seen  at  a  temperature  T  for  all  temperatures  T 
between  Tq  and  Tf.  To  do  this  correctly,  it  must  be  known  how  many  particles  at  any  given 
temperature  T  lie  within  the  sensor  field  of  view  (FOV).  From  Equation  (3.7)  it  is  known  that  all 
the  particles  at  a  given  temperature  T  lie  on  a  sphere  of  radius  r(T)  surrounding  the  flare.  From 
Equation  (3.1),  it  is  known  that  these  exhaust  particles  will  be  evenly  distributed  over  a  spherical 
cap  created,  as  shown  in  Figure  3-4,  by  the  intersection  of  the  conical  exhaust  plume  of  the  flare 
with  a  sphere  of  radius  r(T).  Assume  the  sensor  FOV  is  a  cone  of  half  angle  0^  with  the  burning 
flare  on  the  cone  central  axis  a  distance  s  away  from  the  sensor  (see  Figure  3-1).  The  observation 
angle  0p  is  defined  as  the  angle  between  the  central  axis  of  the  flare  exhaust  plume  and  the  cen¬ 
tral  axis  of  the  sensor  FOV.  When  0p  is  90®,  the  sensor  observes  the  exhaust  plume  broadside; 
when  Op  <  0niax>  sensor  is  located  “inside”  the  flare  exhaust  plume. 

Function  A(r/s,  6p,  6^,  6^^^)  is  defined  as  the  fraction  of  the  total  area  of  the  spherical  cap 
shown  in  Figure  3-4  (located  on  the  surface  of  a  sphere  of  radius  r)  which  lies  inside  the  conical 
sensor  FOV  shown  in  Figure  3-1  (located  a  distance  s  away  from  the  burning  flare).  Obviously, 

A  must  depend  on  dp,  6^,  and  0niax>  t)ut  it  may  not  be  clear  at  first  why  A  depends  only  on  the 
ratio  of  r/s  rather  than  on  r  and  s  as  separate  variables.  Note  that  by  definition,  A  is  a  dimen¬ 
sionless  fraction  which  always  satisfies  the  inequality  0^  A^  1.  Consequently,  distance  can  be 
measured  in  any  convenient  unit  and  there  is  no  reason  not  to  choose  s,  the  distance  between  the 
flare  and  the  sensor,  as  our  unit  of  length.  Thus,  we  can  define  a  new  dimensionless  length 
u  =  r/s  and  know,  because  A,  6p,  6^,  and  6^^^  are  also  dimensionless,  that  A  can  only  depend  on 

0jijax‘ 


Figure  3^4.  How  a  conical  exhaust  plume  intersects  an  imaginary  sphere 
centered  on  the  burning  flare. 
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Function  A  was  calculated  numerically  using  a  Monte  Carlo  computer  program  which 
generated  N  random  points  evenly  distributed  over  a  spherical  cap  a  distance  u  away  from  the 
burning  flare.  (The  size  of  the  spherical  cap,  of  course,  depends  on  both  u  and  ^niax  ) 
of  these  points,  the  computer  program  calculated  an  angle  (see  Figure  3-1)  and  checked 
whether  or  not  that  point  lay  within  the  sensor  FOV  (that  is,  whether  6<}isc  ^  ^s)*  ®y  counting  the 
number  of  points  having  e<jisc  ^  dividing  by  N,  the  program  obtains  an  estimate  for  A. 

The  larger  the  value  of  N,  the  more  accurate  that  estimate.  For  the  results  given  in  this  report, 

1(H. 

A  double  coordinate  transformation  was  used  to  derive  the  formula  for  e^isc-  Th®  central 
axis  of  the  exhaust  plume  is  taken  to  be  the  z  axis  of  the  unprimed  Cartesian  coordinate  system 
shown  in  Figure  3-1.  The  y  axis  of  that  coordinate  system  points  up,  and  the  x  axis  points  out 
of  the  page.  A  single-primed  coordinate  system  is  defined  by  rotating  the  unprimed  system 
through  an  angle  6^  about  its  x  axis  [see  Figures  3-1  and  3-5(a)].  Now  the  z'  axis  points  toward 
the  sensor.  The  second  coordinate  transformation  creates  a  double-primed  coordinate  system  and 
is  a  translation  along  the  z'  axis  by  a  distance  s  [Figure  3-5(b)];  then  the  origin  of  the  double- 
primed  coordinate  system  coincides  with  the  position  of  the  sensor  and  the  negative  z  axis  lies 
along  the  central  axis  of  the  sensor  FOV.  If  x,y,z  are  the  unprimed  coordinates  of  some  point 
lying  on  the  spherical  cap,  and  x'',y'',z''  are  the  coordinates  of  the  same  point  in  the  double- 
primed  coordinate  system,  then 

x''=x, 

y"  =  y  cos0p  +  z  sin0p,  and 

z"  =  -  s  -  y  sinfip  +  z  cos6p  .  (3-11) 

We  can  also  represent  the  position  of  that  point  using  spherical  coordinates  T,6,<f>  in  the 
unprimed  coordinate  system 

X  =  r  sin0  cos<t> 


y  =  r  sin0  sind> 
z  =  r  COS0 


(3.12) 


Similarly,  use  spherical  coordinates  r'',0"’,d)"’  to  represent  the  position  of  the  point  in  the  double- 
primed  coordinate  system.  From  Figure  2-2,  angle  6"  is 


-®disc  =  tan' 


(X")2  +  (y»)2 

z" 


} 


(3.13) 


Hence,  substitute  Equations  (3.12)  into  (3.11),  substitute  the  resulting  equations  into  (3.13),  and 
then  set  ©disc  =  -  6".  This  gives  a  formula  for  ^disc  t®”"®  h  observation  angle),  angles 

e  and  <i>  (the  spherical  angles  in  the  unprimed  coordinate  system),  and  u  (the  dimensionless  dis¬ 
tance  of  the  spherical  cap  from  the  burning  flare). 


(3.14) 
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2  AXIS 

A 


Figure  3-5.  Transformation  from  the  (a)  unprimed  to  the  primed  coordinate  system; 
and  (b)  primed  to  the  double-primed  coordinate  system. 


where 


a  =  sin0  sin<^>  sin0p  -  cos0  cosOp,  and 
b  =  sin^fi  cos^d)  +  (sin0  sind»  cos0p  +  cos©  sin0p)2 

When  the  Monte  Carlo  program  generates  a  random  point,  in  reality  it  chooses  a  pair  of 
random  d,<f>  values  using  the  formula 

e  =  sin-«  j\/l-[l+Q(cos©^ax-l)f  I  (3.15) 

<t>  -  27rQ  . 
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If  Q  is  uniformly  distributed  between  0  and  1,  then  angles  6  and  <f>  will  be  uniformly  distributed 
over  the  unit  sphere  with  0  ^  0  <  e^ax  and  0  <  d>  <  Itt.  Having  generated  a  pair  of  random  e,<t> 
values,  the  program  calculates  from  the  formula  shown  in  Equation  (3.14)  and  decides 
whether  or  not  the  pair  of  d,<f>  values  lie  within  the  sensor  FOV,  Figures  3-6,  3-7,  and  3-8  are 
graphs  of  A  versus  u  for  6^  -  0.1°;  30°,  45°;  and  6p  =  0°,  30°,  60°,  90°. 

In  Figure  3-6,  the  curves  representing  A  for  6p  =  30°,  60°,  and  90°  all  drop  to  zero  at  finite 
values  of  u,  but  the  value  for  6p  =  0°  only  approaches  zero  asymptotically.  This  is  understood  by 
noting  that  for  6p  =  0°  the  sensor  is  always  inside  the  plume  and  will  have  parts  of  the  plume 
within  its  FOV  for  all  values  of  u<l.  In  Figure  3-7,  both  the  curves  for  6p  =  0°  and  for  6p  =  30° 
approach  zero  asymptotically,  and  the  curve  for  dp  =  30°  is  always  less  than  the  curve  for  6p  =  0°. 
Again  this  makes  sense,  because  6p  -  ©nax  "  30°,  so  the  6p  =  30°  curve  lies  on  one  edge  of  the 
plume  as  u  -*  1.  Thus,  it  will  never  reach  zero  and  yet  always  be  smaller  than  the  6p  =  0°  curve, 
which  lies  entirely  inside  the  plume.  In  Figure  3-8,  for  values  of  u  greater  than  0.006  and  less 
than  1,  the  dp  =  30°  and  dp  =  0°  curves  become  identical  (allowing  for  the  randomness  of  the 
Monte  Carlo  estimate).  In  this  case,  dp  <  ~  30°,  so  as  u  —  1  eventu¬ 

ally  both  sensor  locations  will  lie  entirely  inside  the  plume  and  have  the  same  A  values. 


Figure  3-6.  Graph  of  A  versus  u  for  Omax  ”  ■^3°,  6^  =  0.1°,  and  Op  -  0°,  3(F,  60° ,  90° . 


We  define  to  be  the  radiant  energy  per  unit  time  per  unit  wavelength  interval  at  wave¬ 
length  X  per  unit  solid  angle  emitted  by  plume  particles  located  between  r  =  0  and  r  =  Lf.  From 
Equation  (3.5)  and  the  definition  of  A  it  follows  that 


h  =  Go(X,r)A(r/s,6p,0s,e„,ax)  dr  = 


3A^PCi,^(Ro/X) 
47rp  VoRpXS 


T  •A(r/  S,0p,6s,6niax) 

fLfdr-— ^ - 

^0  eC2/AT(r.Ro)_  1 


(3.16) 
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Figure  3-7, 


Figure  3-8. 


Graph  of  A  versus  u  for  6ynax  =  ^0°,  6^  -  0,1^  and  6p  =  (P,  3i 


Graph  of  A  versus  u  for  B^ax  “  Bs  =  O.P  and  Bp  =  (P,  ii 
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Because  A  has  been  evaluated  numerically,  Jx  must  also  be  evaluated  numerically  in  Equa¬ 
tion  (3.16).  It  is  helpful  to  write  this  equation  as 

h 

P 

where  the  variable  of  integration  has  been  changed  from  r  to  T  and  we  have  used  ra(T)  to 
approximate  r(T).  The  formula  for  ldT/dr|  comes  from  Equation  (2.11). 

Equation  (3.17)  is  used  to  calculate  the  data  shown  in  Figures  3-9  through  3-16.  In  these  cal¬ 
culations,  Tq  and  Tf  are  always  set  equal  to  1500®K  and  200®K,  respectively.  Carbon  and  magne¬ 
sium  difluoride  are  the  two  basic  products  of  Mg-Teflon  combustion,  so  p  has  been  set  equal  to 
the  average  of  the  mass  densities  of  bulk  carbon  and  magnesium  difluoride  (2.6  gm/cm^).  The 
value  of  the  heat  capacity  Cp  was  set  to  1.0  X  10^  ergs/gm/®K,  because  both  carbon  (whose  Cp 
values  vary  with  temperature)  and  magnesium  difluoride  have  heat  capacities  of  this  order  of 
magnitude  for  temperatures  between  200°  and  I500°K. 

Figures  3-9  and  3-10  are  graphs  of  Jx/P  for  s  =  100  km;  ~ 

30°,  60°,  and  90°.  [The  “kink”  in  the  graph  of  Figure  3-9  comes  from  the  discontinuous  slope  of 
function,^ defined  in  Equation  (2.18).]  The  ratio  Vol\  300  m/s  (remember  that  \  is  dimen¬ 
sionless)  in  both  Figures  3-9  and  3-10;  their  only  difference  is  that  =  1  /xm  in  Figure  3-9  while 
Rq  =  10  Mm  in  Figure  3-10.  Keeping  in  mind  that  6p  =  90°  is  always  a  broadside  view  of  the 
exhaust  plume,  it  can  be  seen  that  in  each  case  the  amount  of  radiant  energy  Jx  emitted  at  rela¬ 
tively  long  wavelengths  (X  ^  5  Mttt)  increases  as  6p  decreases  to  zero.  This  is  a  purely  geometric 


3Ci,,.^(Ro/X)  _  A(ra(T)/s,  0p,  65,  0|nax) 

- - f  dT - ^ - 

4;rp(Vo/ Ae)RoX5  Tf  IdT/drlfe’^^/xr  _  i) 


Figure  3~9  Graph  of  J)s,IF  versus  X  for  s  =  100  km,  ,  6s-  0,1°, 

VqIA^  -  300  m/s,  Rq-  1  0,437  km. 
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Figure  3-10  Graph  of  J^l P  versus  kfor  s  -  100  km,  Omax  - 
VglAg  =  300  m/s,  Rg  =  10  pm,  and  Lf-  1.92  km. 


phenomenon:  as  6p  gets  smaller,  the  sensor  will  see  more  exhaust  particles  at  greater  distances 
from  the  flare;  being  somewhat  cooler  (but  not  cold),  these  particles  will  present  the  sensor  with 
more  energy  at  longer  wavelengths.  Comparing  Figures  3-9  and  3-10,  it  can  be  seen  that  this 
phenomenon  grows  more  pronounced  as  the  particle  radius  increases.  This  effect  can  also  be 
understood  geometrically.  A  large  particle  cools  more  slowly  than  a  small  one,  thus  making  the 
flare  exhaust  plume  longer  and  bringing  more  cool  particles  into  the  sensor  FOV  at  observation 
angles  6p  <  90°.  Another  way  of  looking  at  this  is  to  note  that  if  all  the  particles  cooled  off 
rapidly  in  the  immediate  vicinity  of  the  burning  flare,  the  plume  would  be  so  short  it  would 
behave  like  a  perfect  point  source,  making  the  Jx/P  curve  look  the  same  for  all  values  of  6p. 

Perhaps  the  best  way  to  analyze  the  increase  of  Jx/P  at  long  wavelengths  is  to  construct  the 
function 


Clearly,  Ex  is  just  the  fractional  increase  of  Jx  at  an  observation  angle  6p  <  90°  compared  to  Jx 
at  an  observation  angle  6p  =  90°  (that  is,  viewed  broadside).  Figures  3-11  through  3-16  graph  Ex 
versus  X  for  6p  =  0°,  30°,  and  60°;  in  each  case  Ex  behaves  as  expected  by  decreasing  as  6p 
increases.  In  all  these  graphs  s  =  100  km,  =  0.1°,  and  Rq  =  10  fim. 

In  Figures  3-11  through  3-13,  v^/Ae  =  100  m/s  and  in  3-14  through  3-16,  \ol\ ~  300  m/s. 
Angle  6njax  ^3°  in  Figures  3-1 1  and  3-14,  30°  in  Figures  3-12  and  3-15,  and  45°  in  Figures  3- 
13  and  3-16.  An  examination  of  Figures  3-11,  3-12,  and  3-13  shows  that  as  0^^  increases.  Ex 


,  =  90°)^  ep=90°) 


(3.18) 
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decreases.  The  easiest  way  to  understand  this  is  to  realize  that  larger  values  of  mean  that 
the  flare  exhaust  plume  fills  up  greater  portions  of  the  spherical  volume  surrounding  the  flare. 
This  gives  the  plume  a  greater  tendency  to  spherical  symmetry,  which  leads  to  a  lesser  depen¬ 
dence  of  Jx  on  0p.  The  same  effect  is  seen  in  Figures  3-14,  3-15,  and  3-16.  Another  point  worth 
mentioning,  after  comparing  Figures  3-11,  3-12,  and  3-13  to  Figures  3-14,  3-15,  and  3-16,  is  that 
as  Vo/Ag  increases,  so  does  Ex-  Larger  values  of  Vo/A^  cause  Lf,  the  exhaust  plume  length,  to 
increase,  and  this  implies  a  greater  geometric  radiance  effect  on  Jx  at  observation  angles 
0p<  90®.  Thus,  increasing  the  value  of  has  the  same  effect  as  increasing  the  value  of  Ro  in 
Figures  3-9  and  3-10;  the  further  the  particles  travel  before  cooling  to  a  given  temperature,  the 
greater  the  dependence  of  Jx  on  6p  at  long  wavelengths. 


5  10  15  20 


X  iftm) 

Figure  3-1}  Graph  of  £\  versus  X  for  s  =  100  km,  Bmax  ^  “  O.P, 

VolAe  =  100  mis,  Rq-  5  pm,  and  Lf  -  0331  km. 


Figure  3-12  Graph  of  Ex  versus  X  for  s  -  100  km,  Bmax  ~  »  B^-  0.1  , 

VolAe  =  100  mis,  Rq-  5  pm,  and  Lf  -  0331  km. 
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Figure  3-13  Graph  of  Ex  versus  k  for  s-  100  km,  Omax  =  45^*  =  0,1 

VqIAc  =  100  mis,  Rq-  5  pm,  and  Lf^  0,331  km. 


Figure  3-14  Graph  of  Ex  versus  X  for  s  =  100  km,  6^ 
VolAff  =  300  m/s,  Rq-  5  pm,  and  Lf-  0.994  km. 


Figure  3-15  Graph  of  £x  versus  k  for  s-  100  km,  0 
Vo/Ag  =  300  m/s,  Ro-  5  pm,  and  Lf-  0.994  km. 


0  5  10 


X  ipm) 


Figure  3-16  Graph  of  £x  versus  kfors-  100  km,  0, 
Vo/Ag  =  300  m/s  Ro  =  5  pm.  and  Lf=  0.994  km. 
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4.  APPLICATION  OF  THE  MODEL  TO  EXPERIMENTAL  DATA 


Figure  4-1  shows  the  basic  experimental  setup  used  for  Bums  6,  7,  and  8  of  the  March  1988 
test  series.  The  flares  were  mounted  at  the  top  of  the  space  chamber  and  aligned  so  that  the  cen¬ 
tral  axis  of  their  exhaust  plumes  pointed  vertically  downward.  Silicon  photometers  (sensitive  to 
wavelengths  between  0.4  and  1.1  iim)  were  mounted  along  the  central  axis  of  the  exhaust  plumes 
at  3  different  stations;  10,  30,  and  91  cm  from  the  flare.  A  2-channel,  Inframetrics  radiometer 
created  a  pair  of  video  images  of  about  the  first  2  m  of  the  flare  exhaust  plume.  The  channel  I 
video  image  came  from  IR  radiation  with  wavelengths  between  2.6  and  3.5  /um,  and  the  chan¬ 
nel  2  video  image  from  IR  radiation  with  wavelengths  between  8  and  12  /zm.  Selected  frames 
were  digitized,  corrected  for  perspective,  and  compared  to  a  calibrated  black-body  source  to  esti¬ 
mate  the  radiant  power  emitted  per  unit  length  of  the  plume  as  a  function  of  distance  down  its 


Figure  4-1.  Experimental  setup  used  in  the  March  1988  test  series. 
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central  axis.  At  the  91-cm  station,  a  helium-cooled,  circular  variable  filter  (CVF)  spectrometer 
viewed  the  exhaust  plume  against  a  liquid-nitrogen-cooled,  cold-wall  background.  The  spectrome¬ 
ter  required  1.2  s  to  scan  the  wavelengths  between  5  and  22.5  fim.  A  typical  bum  lasted  long 
enough  to  provide  between  3  and  5  complete  spectral  scans.  Although  the  formula  for  G(A.,r) 
given  in  Equation  (2.17)  is  for  a  spherical  shell  of  radius  r  (and  thickness  dr)  centered  on  the 
flare,  the  radiometer  estimates  of  radiant  power  are  directly  related  to  G  if,  as  was  the  case  for 
Bums  6,  7,  and  8,  the  plume  spread  angle  e^ax  ‘s  small,  and  r  is  the  distance  down  the  central 
axis  of  the  plume  (see  Figure  2-5).  If  a  function  ^  is  defined  such  that 

»(A„Xb)=  r*‘G(X,r)dX  ,  (“I) 

it  follows  that  channel  1  of  the  radiometer  measured  ^(2.6  /urn,  3.5  (xm)  and  channel  2  measured 
^(8  fim,  12  /Ltm)  as  functions  of  r. 

The  output  of  the  photometer  10  cm  from  the  flare  was  correlated  with  the  output  of  the 
photometer  30  cm  from  the  flare  to  measure  how  fast  the  exhaust  particles  were  traveling  down^ 
the  plume.  The  fluctuating  plume  intensity  seen  by  the  10-cm  photometer  was  labeled  “signal  A”; 
that  seen  by  the  30-cm  photometer  was  labeled  “signal  B.”  In  the  notation  of  our  model,  v^  is 
the  average  velocity  of  the  Mg-Teflon  exhaust  plume  particles.  Fluctuations  in  signal  A  are 
almost  certainly  caused  by  clumps  of  exhaust  particles  traveling  through  the  10-cm  photometer 
FOV,  and  are  usually  followed  by  similar  fluctuations  in  signal  B  after  a  time  lag  r  =  D/v^ 

(where  D  =  0.2  m,  the  distance  between  the  2  photometer  FOVs  along  the  plume  central  aws). 
Thus  finding  the  time  lag  t  which  most  closely  correlates  the  A  and  B  signal  traces  is  equivalent 
to  finding  the  value  of  v^  =  D/t.  Note  that  this  velocity  measurement  automatically  singles  out 
those  particles  which  generate  the  bulk  of  the  radiant  energy  coming  from  the  flare  exhaust 
plume. 

In  Burns  6,  7,  and  8,  0.107-s  matched  sections  of  the  A  and  B  signals  were  sampled  at  17.5- 
fts  intervals.  Thus,  each  0.107  s  section  of  signal  was  turned  into  a  data  set  of  2N  digitized  signal 
values,  N  from  signal  A  and  N  from  signal  B.  Each  data  set  had  a  slightly  different  value  of  N, 
with  6100  <  N  ^  6200  for  every  digitized  data  set  created.  We  define  Aj,  Bj,  for  i  -  0,  I, ... ,  N-1, 
to  be  the  2N  sampled  signal  values  belonging  to  a  given  data  set.  Thus,  Aq,  Bq  would  be  the  first 
A  and  B  signal  samples  of  the  data  set,  Aj,  Bj  would  be  the  A  and  B  signal  samples  17.5  fxs 
later,  and  so  on.  For  each  data  set,  a  cross-covariance  function  Vj  was  defined  such  that 


^  (Ah  -<A»-(Bi-<B» 

1  N-Ij|-1 

^  S  (AHj|-<A».(Bi-<B» 
i=0 


for  j  ^  0 


for  j  <  0 


(4.1a) 
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where 


j  =  -2000,  -1999, ... ,  -1,  0,  1, ... ,  1999,  2000; 


1 


N-l 


i=0 
j  N-l 

<b>=-5;b.  ; 

i=0 

The  cross-correlation  function  Cj  is  defined  in  the  usual  way  as 


Ci=  -7= 

^  \/^o 


where 


N-l 


2  (Ai-<A»2  , 

i=0 

]  N-I 
i=0 


(4.2b) 


When  the  graphs  of  Cj  versus  j  for  all  the  data  sets  were  examined,  about  half  contained 
recognizable  peaks  marking  the  time  lag  which  most  accurately  correlated  the  A  and  B  signals. 
Figure  4-2  shows  the  Cj  versus  j  graph  (from  Bum  8)  which  has  the  most  prominent  peak  found 
in  our  cross-correlation  calculations,  and  Figure  4-3  shows  the  Cj  versus  j  graph  (from  Bum  7) 
which  has  the  least  prominent  peak  that  was  treated  as  a  valid  velocity  measurement.  Table  4-1 
gives  the  average  velocity  measurements  for  Bums  6,  7,  and  8.  Since  roughly  half  the  data  sets 
contain  recognizable  peaks,  this  is  a  good  indication  that  our  model  has  made  a  reasonable 
approximation  of  reality  when  it  assigns  all  the  exhaust  particles  in  the  flare  plume  the  same 
average  velocity  v^,. 

The  output  of  the  10-cm-station  photometer  was  used  to  estimate  the  flare  combustion  rate 
P.  The  area  under  the  photometer  intensity  curve  is  assumed  to  be  proportional  to  the  total 
amount  of  fuel  burned  by  the  flare.  Estimates  were  made  of  the  rectangular  area  that  would  be 
equal  to  the  actual  area  under  each  curve  given  the  constraint  that  the  top  level  of  the  rectangle 
be  equal  in  height  to  the  quasi-steady  intensity  level  characterizing  the  middle  portion  of  each 
burn.  By  dividing  the  length  of  each  rectangle  (representing  a  time  interval  in  seconds)  into  the 
amount  of  fuel  burned  by  each  flare,  estimates  can  be  made  of  the  actual  combustion  rate  P 
characterizing  the  quasi-steady  intensity  levels  of  Burns  6,  7,  and  8  (see  Figure  4-4). 
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Figure  4-2.  Graph  of  Cj  versus  j  from  Burn  8  at  2,9  s  after  ignition. 
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Figure  4-2.  Graph  of  Cj  versus  J  from  Burn  7  at  2.5  s  after  ignition. 


TABLE  4-1 

Average  Velocities  From  Photometer  Correlations 

Over  intervals  of  0.107  s 

Burn  6 
(m/s) 

Bum  7 
(m/s) 

Burn  8 
(m/s) 

Average 

223.5 

238 

273 

Standard  Deviation 

47 

28 

47 

Number  of  Velocity 
Measurements 

6 

4 

9 

A  pyrotechnic  flare  with  a  chemical  composition  similar  to  that  of  Bums  6,  7,  and  8  was 
burned  inside  the  space  chamber  in  such  a  way  that  its  exhaust  plume  passed  through  the  meas¬ 
uring  volume  of  a  Malvern  particle  sizer.  The  particle  sizer  measured  D(4,3)  and  D(3,2)  for  the 
collection  of  particles  in  the  exhaust  plume,  where 


D(4,3)  =  2  •  R4  f(R)dyp  R3  f(R)dR  =  7.4  Mm 

(4.3a) 

D(3,2)  =  2  •  J*®  R3  f(R)dR^  R2  f(R)dR  =  4.9  Mm 

(4.3b) 

Substitute  Equation  (2.2)  into  (4.3a,b)  and  solve  for  a  and  Rj^. 
a  =\/a.[D(4,3)/D(3,2)]  =  0.64 

(4.4a) 

Rla  =  Y  [D(3,2)r/2/[D(4,3)]5/2  =  0.87Mm  . 

(4.4b) 

This  completely  specifies  the  size  distribution  f(R)  defined  in  Equation  (2.2)  and  used  inside  the 
integral  in  Equation  (2.17).  There  is  a  significant  spread  in  this  log-normal  distribution;  the  aver¬ 
age  value  of  R  is  1.1  /xm  while  the  most  probable  value  of  R  is  0.58  fim. 

For  five  different  r  values  (0.1,  0.25,  0.50,  0.75,  and  1  m)  we  calculated  the  average  values  of 
^(2.6  Min;  3.5  Mm)  and  ^(8  Mm,  12  Mm)  using  selected  video  frames  from  the  quasi-steady  portions 
of  Burns  6,  7,  and  8.  These  experimental  values  of  ^  are  presented  in  Tables  4-2,  4-3,  and  4-4. 
Since  the  entire  collection  of  video  frames  was  not  averaged,  when  the  average  of  the  sampled 
frames  is  calculated  it  may,  due  to  the  luck  of  the  draw,  differ  from  the  trae  average  by  a  signifi¬ 
cant  amount.  There  are  well-known  formulas  which  can  be  used  to  calculate  the  standard  devia¬ 
tion  of  this  sort  of  sampling  error,  and  the  percentage  values  of  the  normalized  statistical  error 
given  in  Tables  4-2,  4-3,  and  4-4  are  just  this  standard  deviation  divided  by^.  The  values  of  G/P 
in  these  tables  come  from  the  measured  ^  values  being  divided  by  the  combustion  rate  P  (given 
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BURN  6 


FLARE  WE/GHT==8.94gm;  BURN  RATE  =  P  ^  2.6  gm/s 


BURN  7 


5.4  s 


FLARE  WEIGHT  =  8.16  gm;  BURN  RATE  =  P  =  1 ,6  gm/s 


BURN  8 


FLARE  WEIGHT- 8.90 gm;  BURN  RATE  =  P  -  2.3 gm/s 


Figure  4-4  Graph  of  the  photometer  signal  (lO-cm  station)  versus  time  for 
Burns  6,  7,  and  8. 
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TABLE  4-2 


Average  Radiance  Values  versus  Distance  Down  the  Plume 
for  Burn  6  of  March  1988  Test  Series 

Channel  1  (2.6  ^  X  ^  3.5  /um) 


r  in  m 

^ ,  .  w/sr/m2* 

G/p  in 

kg/s 

^  in  w/sr/m* 

0.1 

2.02X10^0 

47.2 

0.25 

1.61  XIOlO 

37.6 

0.50 

1.11  X  lO^O 

26.0 

0.75 

8.59X109 

20.1 

1.00 

6.67X10® 

15.6 

Channel  2  (8  ^  X  ^  1 2  ^m) 

r  In  m 

^ ,  ,  w/sr/m^* 

G/p  in  .  ^ 

kg/s 

^  In  w/sr/m* 

0.1 

5.88  X  108 

6.11 

0.25 

6.50X108 

6.76 

0.50 

5.75X108 

5.98 

0.75 

4.60  X  108 

4.78 

1.00 

3.77  X  108 

3.92 

•Statistical  error  in  G/p  and  is  ^  4  percent 
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TABLE  4-3 


Average  Radiance  Values  versus  Distance  Down  the  Plume 
for  Burn  7  of  March  1988  Test  Series 


Channel  1  (2.6  <  X  <  3.5  ^m) 


r  in  m 

_  ,  .  w/sr/m2* 

G/P  in  .  , 

kg/s 

^  in  w/sr/m* 

0.1 

2.07X  lOlO 

28.1 

0.25 

1.62X1010 

22.0 

0.50 

1.17X1010 

15.9 

0.75 

9.49  X  109 

12.9 

1.00 

7.73  X  109 

10.5 

Channel  2  (8<X<  12^11) 

r  in  m 

_  ,  w/sr/m^* 

G/p  in  .  , 
kg/s 

^  in  w/sr/m* 

0.1 

7.12  X  IQO 

4.30 

0.25 

7.09X108 

4.28 

0.50 

6.56  X  108 

3.96 

0.75 

4.55  X  108 

2.75 

1.00 

3.99X108 

2.41 

‘Statistical  error  in  G/p  and  ^  is  ^5  percent 
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TABLE  4-4 


Average  Radiance  Values  versus  Distance  Down  the  Plume 
for  Burn  8  of  March  1988  Test  Series 

Channel  1  (2.6  ^  X  ^  3.5  ^m) 


r  in  m 

_  .  .  w/sr/m^* 

G/p  in  .  , 

kg/s 

^  in  w/sr/m* 

0.1 

1.72X1010 

35.7 

0.25 

1.34X  lOlO 

27.7 

0.50 

9.52X109 

19.7 

0.75 

7.15X109 

14.8 

1.00 

5.56  X  109 

11.5 

Channel  2  (8^X^12^m) 

r  in  m 

^ ,  .  w/sr/m^* 

G/p  in 

kg/s 

^  in  w/sr/m* 

0.1 

5.59  X  108 

5.14 

0.25 

5.41  X  10® 

4.98 

0.50 

4.82X108 

4.43 

0.75 

3.43  X  10® 

3.16 

1.00 

2.92X10® 

2.69 

*Statistjcal  error  in  G/p  and  ^  is  <7  percent 
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in  Figure  4-4)  and  the  associated  wavelength  interval  (0.9  /xm  for  channel  1  and  4  for  chan¬ 
nel  1).  The  ^  values  in  Tables  4-2,  4-3,  and  4-4  have  been  plotted  in  Figures  4-5,  4-6,  and  4-7 
with  error  bars  representing  the  un-normalized  standard  deviation  of  the  sampling  error  for  each 
value  of  ^  and  r. 


The  solid  lines  in  Figures  4-5,  4-6,  and  4-7  present  the  best  least-squares  fit  of  the  model 
specified  by  Equations  (2.2),  (2.11),  (2.12),  (2.17),  and  (2.18)  to  the  G/P  values  given  in 
Tables  4-2,  4-3,  and  4-4.  (Finding  the  best  least-squares  fit  to  the  ^  values  themselves  would  have 
taken  an  inordinate  amount  of  computer  time.)  Having  written  a  computer  program  to  calculate 
the  predicted  G/P  values  using  these  5  equations,  it  was  found  that  the  model  values  depend  on 
the  7  free  parameters  Ag,  a,  p,  v^,  Rj^,  Cp  and  T^.  The  2  basic  products  of  Mg-Teflon  combus¬ 
tion  are  carbon  and  magnesium  fluoride,  so  p,  the  mass  density  of  the  exhaust  particles,  is  set 
equal  to  2.6  gm/cm^,  the  average  of  the  mass  densities  of  carbon  and  magnesium  fluoride.  The 
value  of  Vq  comes  from  Table  4-1,  and  the  values  of  a  and  Rja  have  been  calculated  in  Equation 
(4-4a,b).  Thus,  in  the  least-squares  fit,  we  adjust  the  3  remaining  parameters  Ag,  Cp,  and  Tq  to 
fit  the  10  data  points  in  Tables  4-2,  4-3,  and  4-4  to  find  that 


f  I587°K 

for  Bum  6 

To  =i  1561 

7 

(4.5) 

L  1662 

8 

f  1.63X  103  j/kg/°K 

for  Bum  6 

Cp  1.81 

7 

(4.6) 

1  1.20 

8 

r  1.90  (dimensionless) 

for  Bum  6 

Ag  =  ^  2.23 

7 

(4.7) 

11.88 

8 

The  flame  temperature  of  a  chemical  reaction  is  the  predicted  temperature  of  its  products 
given  that  all  the  energy  generated  by  the  reaction  is  absorbed  by  the  products  as  heat.  In  Equa¬ 
tion  (4.5)  the  values  of  Tq,  the  temperature  at  which  the  exhaust  particles  leave  the  flare,  com¬ 
pare  well  to  the  predicted  flame  temperature  (about  2000°K)  of  the  Mg-Teflon  mixture  used  in 
these  3  burns.  On  the  one  hand,  they  are  somewhat  less  than  the  flame  temperature,  which 
makes  sense  because  the  flares  are  not  adiabatically  isolated  thermodynamic  systems;  yet  on  the 
other  hand,  they  are  not  so  much  smaller  as  to  cast  doubt  on  the  model’s  basic  validity.  As  for 
the  values  of  Cp  shown  in  Equation  (4.6),  they  too  seem  reasonable,  given  that  the  Cp  of  pure 
magnesium  difluoride  is  1.21  X  lO^  J/kg/®K  (the  Cp  of  carbon  has  a  pronounced  temperature 
dependence,  yet  at  these  temperatures  its  Cp  is  the  same  order  of  magnitude  as  that  of  magne¬ 
sium  difluoride).  As  for  Ag  in  Equation  (4.7),  it  is  indeed  of  0(1)  as  required  by  our  model. 
Clearly,  the  best  least-squares  fit  of  theory  to  data  does  not  require  the  3  adjustable  parameters 
To,  Cp,  and  Ag  to  take  on  unphysical  values. 


The  error  bars  in  Figures  4-5,  4-6,  and  4-7  represent  the  standard  deviation  of  the  statistical 
sampling  error,  and  it  can  be  seen  that,  by  and  large,  the  channel  1  predicted  radiance  curves  lie 
within  the  error  bars  of  the  experimental  data  for  Bums  6,  7,  and  8.  However  just  the  opposite 
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07853 


Figure  4-5.  Graph  of  average  radiant  intensity  versus 
distance  down  the  plume  for  Bum  6. 
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Figure  4-6.  Graph  of  average  radiant  intensity  versus 
distance  down  the  plume  for  Bum  7. 
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Figure  4-7.  Graph  of  average  radiant  intensity  versus 
distance  down  the  plume  for  Burn  8. 
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can  be  said  about  the  channel  2  predicted  radiance  curves;  they,  as  a  general  rule,  lie  well  outside 
the  error  bars  of  the  experimental  data.  The  best  way  to  understand  this  is  to  remember  that  the 
statistical  sampling  error  is  just  the  random  component  of  the  experimental  error;  the  total  exper¬ 
imental  error  could  be  much  larger  due  to  systematic  inaccuracies  in  the  experimental  apparatus. 
In  fact,  the  general  shape  of  the  channel  2  radiance  data  is  inherently  implausible.  In  Bum  6,  the 
channel  2  radiance  increases  substantially  going  from  r  =  0. 1  to  0.25  m;  and  in  Bums  7  and  8, 
although  the  radiance  decreases,  the  amount  of  the  decrease  seems  anomalously  low.  Further¬ 
more,  although  Table  4-5  shows  that  the  channel  1  radiance  values  are  more  or  less  confirmed  by 
the  photometer  output.  Table  4-6  shows  that  the  CVF  spectrometer  mounted  at  the  91-cm 
station  in  Figure  4-1  is  consistently  in  disagreement  with  the  channel  2  radiance  values.  (The 
channel  2  radiance  is  for  wavelengths  between  8  and  12  nm;  Table  4-6  compares  the  channel  2 
radiance  to  the  CVF  integrated  spectral  intensity  from  8  to  12  ^im.)  If  the  disagreement  between 
the  CVF  spectrometer  and  the  Inframetrics  radiometer  is  taken  to  be  a  measure  of  the  total 
channel  2  error,  the  channel  2  data  are  found  to  be  in  error  by  20  to  30  percent.  If  the  error 
bars  were  to  be  extended  by  this  amount  in  the  channel  2  graphs  of  Figures  4-5,  4-6,  and  4-7,  all 
the  channel  2  experimental  data,  like  that  of  channel  1,  would  match  the  predicted  radiance 
curves.  It  is  thus  tentatively  concluded  that  the  model  radiance  equations  derived  above  ade¬ 
quately  describe,  to  within  experimental  error,  the  radiance  profile  of  Mg-Teflon  flares  burning  in 
a  vacuum. 


TABLE  4-5 

Comparison  of  inframetrics  Radiometer  Channel  1  Readings  to 

91-cm  Station  Photometer  for  March  1988  Test  Series 

Burn 

Number 

Inframetrics 

Channel  1 

Photometer 

1 

0.05  w/cm/sr//Lt 

0.05  w/cm/sr//i 

2 

0.10 

0.08 

3 

0.06 

0.05 

4 

(instrument  failure) 

(instrument  failure) 

5 

0.18 

0.17 

6 

7 

0.13 

0.09 

0.11 

0.08 

8 

0.11 

0.12 

9 

0.16 

0.11 

10 

0.14 

0.12 
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TABLE  4-6 

Comparison  of  Inframetrics  Radiometer  Channel  2  Readings  to 

91 -cm  Station  CVF  Spectrometer  for  March  1988  Test  Series 

Burn 

Inframetrics 

CVF 

Number 

Channel  2 

Spectrometer 

1 

1.1  X  10-^w/cm2/sr/M 

1.1  X  10-4  w/cm2/sr//i 

2 

5.0  X  10-4 

1.9X10-4 

3 

1.3  X  10-4 

0.8X10-4 

4 

{instrument  failure) 

(instrument  failure) 

5 

5.0X10-4 

2.7  X  10-4 

6 

4.4  X  10-4 

3.3  X  10-4 

7 

2.9  X  10-4 

2.3X10-4 

8 

3.3  X  10-4 

2.4X10-4 

9 

3.9  X  10-4 

2.5X10-4 

10 

3.5  X  10-4 

2.0X10-4 
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5.  CONCLUSIONS 


A  model  was  constructed  to  predict  the  IR  radiant  intensity  of  particulate  exhaust  plumes 
generated  by  pyrotechnic  flares  burning  in  a  vacuum.  In  this  model,  it  was  assumed  that  the 
exhaust  plume  particles  come  in  a  variety  of  sizes  described  by  a  log-normal  probability  density 
distribution,  that  the  particles  all  travel  at  the  same  average  velocity,  and  that  the  particles  do 
not  absorb  or  scatter  significant  amounts  of  radiation  —  that  is,  that  the  plume  is  optically  thin. 
Both  the  geometric  irregularity  and  the  lack  of  chemical  purity  of  the  exhaust  plume  particles 
have  been  acknowledged  by  using  a  phenomenological  approximation,  rather  than  Mie  theory,  to 
describe  the  dependence  of  particle  grey-body  emissivity  on  particle  size.  A  single-size  version  of 
this  model  (where  all  the  exhaust  particles  have  the  same  radius  value)  has  been  used  to  predict 
what  the  exhaust  plume  of  a  Mg-Teflon  flare  would  look  like  when  viewed  as  an  approximate 
point  source  by  a  distant  IR  sensor.  The  IR  radiance  data  taken  from  three  separate  Mg-Teflon 
flares  burned  in  the  space  chamber  were  also  compared  to  what  the  multisize  version  of  our 
model  (where  particle  radius  values  obey  a  log-normal  probabilty  density  function)  predicted 
should  occur.  The  single-size  model  gave  a  reasonable  description  of  a  Mg-Teflon  flare  viewed  as 
a  distant  source  of  IR  radiation,  and  the  multisize  model  predictions  matched  the  space-chamber 
data  to  within  estimated  experimental  error.  Clearly,  the  approximations  used  to  construct  the 
model  are  a  reasonable  description  of  reality,  and  the  model  itself  is  a  useful  tool  for  understand¬ 
ing  the  exhaust  plumes  of  pyrotechnic  flares  burning  in  a  vacuum. 
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